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Abstract 

Self-incompatibility (SI) of tlie Brassicaceae family can be overcome by CO2 gas treatment. This metliod has been used 
for decades as an effective means to obtain a large amount of inbred seeds which can then be used for hybrid seed 
production; however, the molecular mechanism by which CO2 alters the SI pathway has not been elucidated. In this study, 
to obtain new insights into the mechanism of C02-induced SI breakdown, the focus was on two inbred lines of Brassica 
rapa (syn. campestris) with different CO2 sensitivity. Physiological examination using X-ray microanalysis suggested that 
SI breakdown in the C02-sensitive line was accompanied by a significant accumulation of calcium at the pollen-stigma 
interface. Pre-treatment of pollen or pistil with CO2 gas before pollination showed no effect on the SI reaction, suggesting 
that some physiological process after pollination is necessary for SI to be overcome. Genetic analyses using progeny 
of a C02-sensitivexC02-insensitive cross suggested that CO2 sensitivity is a semi-dominant trait in these lines. Analysis 
of F2 progeny suggested that CO2 sensitivity could be a quantitative trait, which is controlled by more than one gene. 
Quantitative trait locus (QTL) analyses identified two major loci, BrSIOI and BrSI02, which work additively in overcoming 
SI during CO2 treatment. No QTL was detected at the loci previously shown to affect SI stability, suggesting that CO2 sensi- 
tivity is determined by novel genes. The QTL data presented here should be useful for determining the responsible genes, 
and for the marker-assisted selection of desirable parental lines with stable but C02-sensitive SI in F^ hybrid breeding. 

Keywords: Brassica rapa, calcium, CO2, hybrid, QTL, self-incompatibility. 



Introduction 

Self-incompatibility (SI) is a widespread genetic system in 
many flowering plants which serves to prevent self-fertiliza- 
tion and maintain genetic diversity. It is based on self/non-self 
pollen-pistil recognition interactions followed by inhibition 
of self-pollen hydration, germination, or pollen tube growth. 

In the Brassicaceae, SI is sporophytically controlled by a 
multialleHc locus termed the S locus (Bateman, 1955). Male 
and female determinants have been identified as SPll/SCR 
(5'-locus protein ll/5'-locus cysteine-rich) (Schopfer et ciL, 
1999; Takayama et al, 2000) and SRK {S receptor kinase) 



(Takasaki et al, 2000), respectively. When a compatible pol- 
len grain lands on the stigma, it swells and a pollen tube is 
allowed to grow, whereas when self-pollen attaches to the 
stigma, SPll/SCR binds specifically to the extracellular 
domain of SRK of the same S-haplotype (Takayama et al. , 
2001), which triggers an SI signalHng pathway to reject self- 
pollen. Another stigmatically expressed gene located at the 
5" locus, S locus glycoprotein (SLG) (Nasrallah et al., 1987; 
Takayama et al., 1987), has been shown to enhance the rec- 
ognition process between self-pollen and stigma (Takasaki 



Abbreviations: InDel, insertion/deletion; LG, linkage group; LOD, log of odds; QTL, quantitative trait locus; RFLR restriction fragment length polymorphism; 

RLSICO2, reaction level of SI to CO2; SI, self-incompatibility; SSR, simple sequence repeat. 
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et al., 2000); however, this function of SLG remains con- 
troversial (Silva et al., 2001). A recent study suggested that 
the plants in the Brassicaceae genus Levenworthia use par- 
alogous SRK and SPlllSCR genes, Lal2 (Leavenworthia ala- 
hamica SRK-related 2) and SCRI(SCR-like), for self/non-self 
recognition in SI, but the function of their orthologues in 
other Brassicaceae genera also remains unknown (Chantha 
et al, 2013). 

Many studies focusing on the downstream components 
involved in this type of SI signalling pathway have been per- 
formed and, thus far, two components have been identified as 
positive effectors. ARCl (arm repeat containing 1) was iden- 
tified by a yeast two-hybrid screen using the kinase domain of 
SRK as the bait (Gu et al., 1998; Stone et al., 1999). ARCl is 
a U-box protein with E3 ubiquitin ligase activity (Stone et al. , 
2003), and has been shown to interact with ExoVOAl, a puta- 
tive component of the exocyst complex required for com- 
patible pollination (Samuel et al., 2009). MLPK (M-locus 
protein kinase) was identified by positional cloning as the 
gene responsible for the self-compatible mutation of Brassica 
rapa var. Yellow sarson (Murase et al., 2004). MLPK is a 
membrane-anchored cytoplasmic protein kinase and inter- 
acts directly with SRK to transduce SI signalling (Kakita 
et al, 2007). However, the importance of these two compo- 
nents in the Brassicaceae SI mechanism remains controversial 
(Kitashiba et a/., 2011; Indriolo et al, 2012). 

In the Brassicaceae, it has been known that SI can be overcome 
under some physiological and environmental conditions such 
as plant age (Ockendon, 1978; Horisaki and Niikura, 2008), 
stigmatic chemical treatments (e.g. ether, KOH, and NaCl) 
(Tatebe, 1968; Tao and Yong, 1986; Monterio and Gabelman, 
1988), and high temperature (Matsubara, 1980; Okazaki and 
Hinata, 1987). CO2 gas (3-5%) treatment (Nakanishi et al, 
1969) is the most effective way to overcome SI. Today, most 
cultivated lines of crucifer vegetables, such as cabbage, broc- 
coli, Chinese cabbage, and radish are Fj hybrids whose seeds 
are produced by mix-planting two self-incompatible inbred 
parental lines. In this economical Fj hybrid breeding system, 
CO2 gas treatment has been used to suppress SI and allow self- 
fertilization, thereby providing large-scale seed propagation of 
parental lines. This method has been used all over the world for 
many years, but the molecular mechanism leading to SI break- 
down by CO2 gas treatment is entirely unknown. 

Previous studies suggested that not all lines respond equally 
to CO2, and there are variations in SI response to CO2 (CO2 
sensitivity) in the Brassicaceae (Nakanishi and Hinata, 1973; 
Niikura and Matsuura, 2000). Preliminary genetic analysis 
using lines with different CO2 sensitivity in radish (Raphanus 
sativus) suggested that high CO2 sensitivity was controlled 
by a recessive gene independent of the ^-locus (Niikura 
and Matsuura, 2000). In another genetic analysis using 
C02-sensitive and C02-insensitive lines of Chinese cabbage 
{B. rapa), high CO2 sensitivity was suggested to be controlled 
by a dominant gene (Hyun et al. , 2007); however, no responsi- 
ble genes have been identified from these studies so far. 

In this study, new inbred lines of B. rapa with different 
CO2 sensitivity, a C02-sensitive line (HA-11621) and a CO2- 
insensitive line (HA-11623), were selected and analysed. 



X-ray microanalysis suggested that SI breakdown in the 
C02-sensitive line was accompanied by significant calcium 
accumulation at the pollen-stigma interface. Independent 
pre -treatment of pollen or pistil with CO2 gas before pollina- 
tion showed no effect on the SI reaction, suggesting that some 
physiological process that occurs after pollination is neces- 
sary for SI to be overcome. Genetic analyses using Fi and F2 
progeny of a C02-sensitivexC02-insensitive cross suggested 
that CO2 sensitivity is a semi-dominant and quantitative trait. 
Furthermore, quantitative trait locus (QTL) analyses identi- 
fied two major responsible loci, BrSIOl and BrSI02, which 
function additively in overcoming SI during CO2 treatment. 

Materials and methods 

Plant materials 

Two inbred lines of B. rapa (2/2=20), a C02-sensitive line (HA- 
11621) and a C02-insensitive line (HA-11623), were established at 
Tohoku Seed Co., Ltd, and grown in the greenhouse with 16h light 
and 8 h dark conditions at 20 °C. Both lines show stable SI under 
normal (open-air) condition but have different CO2 sensitivity; SI in 
HA-1 1621 breaks down following treatment with 4.5% CO2 whereas 
SI in HA-1 1623 is unaffected. HA-1 1621 and HA-1 1623 are recipro- 
cally compatible, and their Fj progeny were obtained under normal 
conditions. Buds (1-2 d before flowering) from a randomly chosen 
F] were used for F2 production. Young petals and stamens were 
removed from the bud, and the immature pistil was pollinated with 
pollen grains from mature flowers of the same plant (bud pollina- 
tion). Pollinated pistil was then covered with a paper bag for 3 d and 
seeds from the pistil were harvested. More than 20 pistils were pol- 
linated, and harvested seeds were used as the F2 population. A total 
of 1 10 F2 plants were used for phenotypic and genetic analysis. 

Cryo-scanning electron microscopy and energy-dispersive 
X-ray analysis 

Flowers were self- or cross-pollinated and incubated for 1.5 h with or 
without 4.5% CO2 gas. These pollinated and non-pollinated pistils 
were submersed in liquid nitrogen slush and frozen under vacuum. 
While under vacuum, the sample was transferred to the microscope 
cryo stage (ALTO 1000, Gatan), and then the stage temperature was 
increased to -95 °C to remove frost that had settled on top of the 
specimen as a result of condensation. When all surface frost had been 
removed by sublimation, as verified by electron microscopy, the tem- 
perature was reduced to -140 °C. Imaging was performed using an 
ETD (Everhart-Thornley detector) by Quant 250 scanning electron 
microscopy (FEI). The chamber pressure was 30 Pa and the acceler- 
ating voltage was 1 5 kV. EDX (energy-dispersive X-ray spectroscopy) 
analysis of the element assay was performed on selected papilla cells 
using INCA X-ray analysis software (Oxford Instruments, http:// 
www.oxinst.com/Pages/home.aspx, last accessed 14 December 2013), 
with the detector's processing time set at 2. X-ray data were coUected 
with 4.5 nA probe current for 2min. Each 2-3 pistils were used in one 
experiment and three individual experimental sets were performed. 

Evaluation of reaction level of SI to CO2 (RLSICO2) 

Three to five flowers were cut at the peduncle and stood on a 1% 
(w/v) solid agar plate. Flowers were self-pollinated, placed into a 
CO2 incubator, and treated with 4.5% CO2 for 4h at 23 °C. After 1 
d at room temperature, pistils were fixed in ethanol: acetic acid (3:1) 
overnight, softened in 1 N NaOH at 60 °C for 2h, then stained with 
0.01% (w/v) decolorized aniline blue in 2% K3PO4 for 6h. Pollen tube 
behaviour was observed under a fluorescent microscope (Axiophot 
2, Zeiss). CO2 sensitivity was measured using the RLSICO2 index. 
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RLSICO2 was classified into five categories, based on the number of 
pollen tubes penetrating into the stigma: 1, 0 pollen tubes; 2, 1-5 pol- 
len tubes; 3, 6-15 pollen tubes; 4, 16-30 pollen tubes; and 5, >30 pol- 
len tubes. Three replicates were performed on each plant on different 
days. Non-C02-treated self-pollinated flowers were used as controls. 
In aU cases, no pollen tubes penetrated into the control stigmas. 

Genotyping of S-haplotypes 

S'-haplo types of B. rapa were identified using primers PS 5 
(5'-ATGAAAGGCGTAAGAAAAACCTA-30 and PS15 (5'-CCG 
TGTTTTATTTTAAGAGAAAGAGCT-3') (Nishio et a!., 1996) 
to amplify a fragment of the SLG gene. PCR-RFLP (restriction 
fragment length polymorphism) was used to distinguish the two 
S-haplotypes based on differential digest with the restriction enzyme 
Kpnl (TaKaRa, Japan). Digested DNA was electrophoresed on a 
1.5% agarose gel. 

Molecular markers and detection of DNA polymorphism 

To screen for markers that show polymorphism between B. rapa 
lines, primers specific for simple sequence repeat (SSR) markers from 
different sources [UK, prefixes Ra, Na, Ol, and ENA (Lowe et al. , 
2004; http://brassica.bbsrc.ac.uk); Japan, prefixes BRMS, KBr, and 
EST (Suwabe et al, 2002, 2004, 2006; http://vegmarks.nivot.affrc. 
go.jp, NIVTS); China, prefix sau_um (Ge et al, 2011); and Korea, 
prefix AMCP (Ramchiary et al. , 201 1)] were used. SSR, RFLP, and 
insertion/deletion (InDel) markers (prefixes XT and Bra) were also 
designed based on the Brassica database (BRAD) (http://brassicadb. 
org/brad/, last accessed 14 December 2013) (Supplementary Table 
SI available at JXB online). 

Total genomic DNA was extracted from young leaves of two 
parental lines and F2 progeny using the cetyltrimethylammonium 
bromide (CTAB) method (Murray and Thompson, 1980). DNA 
polymorphism analysis with SSRs was carried out using PGR with 
fluorescent dyes, performed according to Suwabe et al. (2008) with 
some modifications. The M13 (-21) universal primer sequence 
(18 bp) was fused to the 5' end of the original forward primer, and 
the M13 (-21) universal primer was labelled with 6-FAM, NED, 
VIC, or PET fluorescent dye (Applied Biosystems, CA, USA). PCRs 
were performed in a 10 |il reaction volume containing lOng of tem- 
plate DNA, 4.7 |xM of labelled M13 (-21) universal primer and 
reverse primer, 0.3 |xM of forward primer, 1 x PCR buffer, 1 x dNTP, 
ix MgCl2, and 0.5 U of rTaq (TOYOBO, Japan). Conditions for 
PCR were as follows: initial denaturation was carried out at 94 °C 
for 3min foOowed by 37 cycles at 94 °C for 30 s, 55 °C (slope of 
0.5 °C s-') for 30 s, 72 °C (slope of 0.5 °C s"') for 30 s, and a final 
extension at 72 °C for 4min. A 1 \i\ aliquot of 50-fold diluted PCR 
product was added to 8.9 |il of Hi-Di™ Formamide and 0.1 |il of 
GeneScanTM 600 LIZ™ Size Standard (Applied Biosystems, USA) 
and applied to an ABI 3730 DNA Analyzer (Applied Biosystems). 
Data were analysed using ABI GeneMapper* software. 

For polymorphism analysis with RFLP and InDel markers, PCR 
was carried out in a 10 |il reaction volume with 5 pmol of forward 
and reverse primers instead of fluorescent dyes. For RFLP markers, 
amplified fragments were digested using restriction enzymes for 1 h. 
Fragments of digested DNA were separated on a 2-4% agarose gel. 

Linl<age map construction and QTL analysis 

A genetic map was constructed using JoinMap® version 4 (Van 
Ooijen, 2006) utilizing the double pseudo-testcross strategy with a 
logio of odds (LOD) threshold of 6.0 for linkage group identifica- 
tion. The best marker order was calculated with the regression map- 
ping algorithm, and marker order was retained from the first round 
only. Map distance units in centiMorgans (cM) were converted from 
recombination frequencies using the Kosambi mapping function 
(Kosambi, 1944). Interval mapping (IM) was performed to identify 
putative QTLs using the established linkage map and the observed 



phenotypic traits. This method was run using MapQTL® version 6 
(Van Ooijen, 2009). With this software, a P < 0.05 LOD score signifi- 
cance threshold was calculated by creating a group-wide distribution 
of the data based on a 1000 permutation test. LOD peaks were used 
to estimate the position of QTLs on the map. 

Statistical analysis 

Box plots were prepared by Ekuseru-Toukei 2012 software (Social 
Survey Research Information Co., Ltd, Japan) to compare the phe- 
notypic difference, as this plot type gives a good sense of environ- 
mental data distribution (Upton and Cook, 1996). Kruskal-Walhs 
analysis of variance (ANOVA) by ranks was used between paired 
comparisons of markers to examine marker association. 

Results 

Phenotypic analysis of B. rapa lines in response to 
CO2 treatment 

Two inbred lines of B. rapa with different CO2 sensitivity, a 
C02-sensitive line (HA- 11 62 1 ) and a C02-insensitive line (HA- 
11623), were used in this study. Flowers were self-pollinated 
by hand pollination and incubated in a CO2 incubator (4.5% 
CO2) for 4h. Both lines were self-incompatible under normal 
conditions (control) (Fig. 1 A, B), whereas they showed signif- 
icantly different responses to CO2 gas treatment (Fig. IC, D). 
Specifically, in the C02-sensitive line, many pollen tubes were 
seen to penetrate into papilla cells after CO2 treatment. This 
pollination test confirmed that the C02-sensitive line had 
high sensitivity to CO2 while the C02-insensitive line hardly 
responded to 4.5% CO2. Cross-pollination was performed as 
a positive control (Fig. IE, F). 

Physiological changes in papilla cells after CO2 
treatment 

Previous work using X-ray microanalysis has revealed the 
accumulation of calcium at the stigmatic surface following 
compatible cross-pollination in Brassica oleracea (EUeman 
and Dickinson, 1999). X-ray mapping strongly indicated 
that a high concentration of calcium was localized at the 
points where the pollen grain made contact with the surface 
of the stigmatic papilla cell. The calcium accumulation was 
also observed in B. rapa, especially in compatible pollination 
(Iwano et al, 1999). To examine the physiological effect of 
CO2 on pollination reactions, this calcium accumulation was 
examined by using cryo-scanning electron microscopy fitted 
with an X-ray microanalysis system. When the C02-sensitive 
line was cross-pollinated with the C02-insensitive line, pollen 
grain hydrated and the pollen tube penetrated into the papilla 
cell, but when C02-sensitive and C02-insensitive lines were 
self-pollinated under normal conditions (without CO2 treat- 
ment), few pollen grains hydrated, and no pollen tube ger- 
mination was observed in either line (Fig. 2A, upper panel). 
After CO2 treatment, the cross-pollinated pollen grains did 
not show a significant difference, neither did the self-polli- 
nated C02-insensitive line. However, the self-pollinated CO2- 
sensitive line showed obvious changes: pollen hydration and 
germination were observed under CO2 treatment (Fig. 2A, 
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Fig. 1. Phenotype of inbred Brassica rapa lines used in tinis study. (A, B) Pollen tube behaviour after self-pollination of COj-sensitive 
(HA-11621) and C02-insensitive (HA-1 1623) lines under normal conditions (without COg treatment). No penetrated or elongated pollen 
tubes are observed in either line. (C, D) Pollen tube behaviour after self-pollination of COj-sensitive and C02-insensitive lines under 4.5% 
CO2 gas treatment. Pollen tubes could penetrate into the stigma and elongate through the style only in the C02-sensitive line. (E, F) 
Cross-pollination as a positive control; the arrow shows pollen tubes which have penetrated. Pt, pollen tubes. Bar=1000 nm. 



lower panel). The emission of elements (Ka) in the tip of the 
papilla cells was then analysed after selecting the pollinated 
papilla cells which faced nearly the same direction in relation 
to the X-ray detector. Emissions of P-Ka, S-Ka, K-Ka, and 
Ca-Ka were detected together with C and O, which are con- 
stituent elements of biological materials. The detected emis- 
sion of Al-Ka originates from the stub that held the samples 
(Fig. 2B). There was no significant difference in the elemental 
emissions between C02-sensitive and C02-insensitive lines 
before pollination. Ca-^ accumulation was observed after 
cross-polHnation as previously reported (Iwano et al., 1999), 
and Ca-Ka emission was increased with CO2 treatment. 
After self-pollination, Ca-Ka emission was slightly increased 
in the C02-sensitive line in the normal CO2 condition. This 
Ca-Ka increase became more pronounced (~6-fold) in the 
C02-sensitive line when it was self-pollinated in the high CO2 
condition. In the C02-insensitive line, no significant Ca-Ka 
increase was observed after self-pollination in normal or high 
CO2 conditions (Fig. 2B). Although the biological signifi- 
cance of the accumulation of calcium at the pollen-stigma 
interface is not clear, it has been suggested that calcium plays 
some role in the successful development of the pollen tube 
tip into the region of expanded stigmatic wall (Elleman and 
Dickinson, 1999). The present results suggest that, in the 
C02-sensitive line, high CO2 activates a compatible pollina- 
tion pathway or disturbs an SI signalling pathway leading to 
Ca^"^ accumulation at the pollen-stigma interface. 

The efficiency of CO2 treatment 

A previous study suggested that the effect of CO2 on SI 
breakdown depends on the timing of treatment (Nakanishi 



and Hinata, 1973). In the C02-sensitive Hne, when self-polH- 
nated flowers were immediately treated with 4.5% CO2 for 4h, 
SI could be overcome, and typically >10 pollen tubes pene- 
trated into the stigma (Fig. 3A). When CO2 treatment started 
3h after self-pollination, SI could still be overcome (Fig. 3B). 
However, when CO2 treatment started 6h after self-pollina- 
tion, the number of penetrating pollen tubes was decreased 
(Fig. 3C). These results indicate that self-pollen inhibition in 
SI is biostatic, as previously suggested (Sarker et al., 1988), 
and can be reversed at least at 3 h after pollination. However, 
at 6h after pollination, SI inhibition enters an irreversible 
phase that cannot be overcome by CO2 treatment. 

Next, in order to narrow down the stage of SI affected by 
CO2 treatment, experiments with high C02-pre-treated pollen 
or pistil from the C02-sensitive line were performed (Fig. 4). 
SI could not be overcome by pre-treatment of either tissue, 
and no pollen tube penetration could be observed even when 
both pollen and pistil were treated separately prior to pollina- 
tion (Fig. 4). SI could be overcome only when the self-polH- 
nated pistil was treated with high CO2. These results suggest 
that some post-pollination physiological process is affected 
by high CO2, in the process of SI breakdown. 

S-allele cliaracterization and plienotype of CO2 
sensitivity in F, and F2 

The S-haplotypes of the two parental inbred lines were first 
determined by ampHfying their SLG genes (Nishio et al, 
1996). The sequence data suggested that the 5-haplotypes 
of the C02-sensitive and C02-insensitive lines were 5*55555 
and 546546, respectively. To dissect genetically the gene(s) 
that determines sensitivity to CO2 treatment, six Fj plants 
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Fig. 2. Electron micrographs and X-ray microanalysis of S. rapa. (A) Cryo-scanning electron micrograplns of pollinated papilla cells 
were taken 1 .5h after pollination. Representative examples of the cross-pollinated (left column) and self-pollinated (middle column) 
COg-sensitive line, and the self-pollinated (right column) C02-insensitive line are shown. Without CO2 treatment (upper panels), only 
cross-pollen is accepted, and self-pollen grains maintain a spheroid shape without swelling in both lines. With 4.5% CO2 gas treatment 
(lower panels), pollen grains swell and germinate in the C02-sensitive line (lower middle) but not in the C02-insensitive line (lower right). 
Bar=25 \im. (B) Representative examples of energy-dispersive X-ray spectra of non-pollinated and pollinated papilla cell surfaces. 
Scanning positions for X-ray analyses are indicated by asterisks in (A). The emissions of Al-Ka, P-Ka, S-Ka, K-Ka, and Ca-Ka were 
detected at the papilla cell surfaces, and the intensity of Ca emission was increased after cross-pollination (left column). The increase 
of Ca emission was also observed after self-pollination with 4.5% CO2 gas treatment in the C02-sensitive line (middle column) but not 
in the C02-insensitive line (right column). These spectrum patterns are reproducible in three individual experiment sets. The emission of 
Al-Ka is mostly derived from the stub that held the samples. 
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Fig. 3. The efficiency of COg treatment of the COg-sensitive line. (A) Self-pollinated stigmas were treated with 4.5% COg immediately 
after self-pollination, and pollen tube penetration and elongation were observed. (B) Self-pollinated stigmas were treated with 4.5% COg 
3h after self-pollination. The pollen tube can still penetrate into papilla cells. (C) When CO2 treatment began 6h after self-pollination, the 
number of penetrating pollen tubes was decreased. Pt, pollen tubes. Bar=1 000 nm. 
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Fig. 4. Pollination assay using C02-treated non-pollinated flowers of the C02-sensitive line as either COj-treated pistil or pollen. Pollen 
tube behaviour is observed after a 4 h CO2 treatment immediately after self-pollination. (A) A pistil from a non-pollinated flower treated 
with CO2 for 4h was pollinated with pollen from the same plant which was not treated with CO2. (B) A pistil from a non-treated non- 
pollinated flower was pollinated with pollen from the same plant which was treated with CO2 for 4h. (C) A pistil from a C02-treated 
non-pollinated flower (C02-treated) was pollinated with pollen from the same plant which was treated with CO2 for 4h. (D) Self-pollinated 
flower treated with CO2 after pollination, positive control. Pt, pollen tubes. Bar=1000 |.im. 



(546555) were produced by crossing C02-sensitive and CO2- 
insensitive lines. These Fj plants exhibited an intermediate 
CO2 sensitivity phenotype where self-pollen tubes penetrated 
into the stigma under high CO2 treatment but there were 
fewer penetrating pollen tubes than observed in a self-polH- 
nation of the C02-sensitive parent. An F2 population of 110 
individuals derived from a bud-pollinated Fi plant was made 
and used for further genetic analyses of the CO2 sensitivity. 
F2 individuals were genotyped using PCR-RFLP to distin- 
guish SLG alleles (Supplementary Fig. SI at JXB onHne). 
555- and 546-haplotypes were segregated in the F2 popula- 
tion according to MendeHan transmission (Supplementary 
Table S2). Pollen tube behaviour after CO2 treatment varied 
among individuals and, in order to quantify the strength of 
CO2 sensitivity, the modified RLSICO2 (reaction level of SI 
to CO2) index was employed (Niikura and Matsuura, 2000), 
which calculates CO2 sensitivity based on the number of pen- 
etrating pollen tubes after self-pollination under high CO2 



conditions (see the Materials and methods). The RLSICO2 
of 110 F2 individuals is presented in Supplementary Fig. 
S2, and the summarized box-plot data are shown in Fig. 5, 
together with the RLSICO2 of Fi and the parental inbred 
fines. Fi had an RLSICO2 score intermediate to the two 
parental lines, suggesting tfiat the high CO2 sensitivity is a 
semi-dominant (incompletely dominant) trait in these inbred 
fines. Furthermore, the RLSICO2 of F2 individuals was con- 
tinuously distributed and did not follow a simple one-locus 
biafielic Mendefian distribution (Supplementary Fig. S2). 
These results suggest that CO2 sensitivity in the inbred lines 
used here could be a quantitative trait which is controlled by 
more than one gene. 

Relationship between S-alleles and CO2 sensitivity 

To investigate whether CO2 sensitivity is related to 
5-haplotypes, the 110 F2 individuals were grouped into three 
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genotypes (S55S55, 84^84^, and 84^85^)- The RLSICO2 of each 
group is shown in Fig. 5. In the three F2 groups, RLSICO2 
scores were distributed from 1 to 5, and interquartile ranges 
overlapped, indicating that CO2 sensitivity is not linked to the 
S'-locus in these two lines. 



Reproductive tissue controiiing CO2 sensitivity 

From the F2 population, two 84^84,, homozygotes with dif- 
ferent RLSICO2 were selected: F2-I6, a C02-insensitive 
line (RLSIC02=1 ±0); and F2-26, a CO^-sensitive line 
(RLSICO2=4.15±0.53). These two lines were used to 
examine the reproductive tissue controlling CO2 sensitiv- 
ity. Reciprocal crosses were performed between these two SI 
lines with or without high CO2 gas treatment. Because all 
crosses under normal conditions (without CO2 treatment) 
were incompatible, only data from crosses performed in the 
high CO2 condition are shown in Fig. 6. The cross between 
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I 

COTinsen 



■^46'^55 
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Fig. 5. Box plots of CO2 sensitivity phenotypes. Data slnow the 
distribution of RLSICO2 witli 25tln, 50tli, and 75th percentiles 
(horizontal bars), interquartile ranges (columns), and 1.5 
interquartile ranges (error bars) of RLSICO2 from six C02-sensitive 
(S55S55) and six C02-insensitive (S46S46), individuals, six 
individuals (S46S55), and 110 F2 individuals (22 S46S46, 64 S46S55, 
and 24 S55S55). a indicates a significant difference (P < 0.01) 
between C02-sensitive and C02-insensitive lines. 



C02-sensitive F2-26 pistil and C02-insensitive F2-I6 pollen 
was CO2 sensitive, showing many penetrating pollen tubes 
under high CO2. On the other hand, the cross between CO2- 
insensitive F2-I6 pistil and C02-sensitive F2-26 pollen was 
CO2 insensitive, showing no penetrating pollen tubes even 
under high CO2. These results suggest that CO2 sensitivity 
is controlled by genes expressed in the female tissue (pistil). 

Mar/cer anaiysis and construction of a linkage map 

In order to map QTLs that determine CO2 sensitivity, a link- 
age map was constructed for this F2 population. A total of 
911 different genetic markers were examined in the two paren- 
tal lines. To clarify the relationship between previously iden- 
tified Sl-related genes and CO2 sensitivity, 8LG (an S-locus 
marker), MLPK, and ARC! were also selected. Though a 
very low level of polymorphism (14.7%) was detected for all 
types of markers, 123 polymorphic markers were selected, 
which include 113 SSRs, five RFLPs, and five InDel mark- 
ers. These 123 markers were used for linkage mapping, and 
generated 10 Hnkage groups (AOl-AlO) at a LOD threshold 
value of 6.0 (Fig. 7). The total length of the map was 947.5 
cM, and the length of the linkage groups ranged from 63.2 
cM (AlO) to 168.4 cM (A03). The distance between markers 
varied from 0 to 29.3 cM, with an average interval of 7.7 
cM. 8LG, MLPK, and ARCl were mapped to A07, A03, and 
A04, respectively, which is consistent with previous reports 
(Ajisaka et al, 2001; Hatakeyama et al, 2010). 

QTL anaiysis and association of markers witii higii CO2 
sensitivity 

Using the constructed linkage map, QTLs responsible for high 
CO2 sensitivity were analysed. Three QTLs were identified on 
linkage groups A03 and AOS based on a LOD threshold of 
3.40 (1000 permutation test, P < 0.05) (Fig. 8, Table 1). These 
QTLs were tentatively named Brassica rapa 81 Overcome 
( Br8I0 ) 1-3, and these results further supported the predic- 
tion that CO2 sensitivity of SI is controlled by a polygenic 
system. Br8I01 on AOS and Br8I02 on A03 are two major 
QTLs that explained 19.3% and 19.0% of phenotypic varia- 
tion, respectively. Br8I03, located near Br8I02, accounted 
for 14.5% of the variance (Table 1). 
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Fig. 6. Reciprocal crosses with CO2 treatment between two S46 homozygous individuals from the F2 population with different RLSICO2. 
(A) C02-sensitive F2 self-pollination. (B) C02-insensitive F2 self-pollination. (C) A C02-sensitive F2 pistil pollinated with pollen from a CO2- 
insensitive F2. (D) A C02-insensitive F2 pistil pollinated with pollen from a C02-sensitive F2. Pt, pollen tubes. Bar=1000 |im. n=3. 
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Fig. 7. A linkage map of selected DNA markers from the population. IVlap distances are shown to the left of vertical lines of the 
linkage group in cM, and marker names are shown to the right. 



To examine the significance of these QTLs, F2 progeny were 
classified into groups based on the genotypes of the hnkage 
markers nearest these three newly identified loci, and the 
relationship of the loci to RLSICO2 in individual plants was 
analysed using Kruskal-Wallis ANOVA by ranks (Table 2). 
Alleles from C02-sensitive (HA- 11 621) and C02-insensitive 
(HA- 11 623) lines are presented as S and /, respectively. 
Almost all classifications using the closest linkage markers 
showed a higher RLSICO2 index in the SS group with signifi- 
cance dt P < 0.01, except marker BRMS-114, which showed 
significance at P < 0.05. 

Marker association was further examined with combina- 
tions of BrSIOI and BrSI02 since BrSI03 is a minor QTL 
closely linked to BrSIOl, making it difficult to identify it 
as an independent QTL. F2 progeny were classified into 
nine groups based on the genotypes of their closest link- 
age markers (Table 3). According to this classification, for 
example, the two above analysed Si,(,S^(, homozygous lines 
with different RLSICO2, F2-I6 (C02-insensitive line) and 
F2-26 (C02-sensitive line), were classified into group 8 and 
group 1, respectively. When groups with the same BrSIOl 
genotype were compared, the BrSI02 SS group showed a 
higher RLSICO2 index compared with the BrSIOl II group. 
Likewise, when the groups with the same BrSI02 genotype 
were compared, the BrSIOI SS group showed a higher 



RLSICO2 compared with the BrSIOI //group. Although the 
numbers of F2 individuals in each group were rather low, sig- 
nificance (P < 0.05) was detected between groups 1 and 6, 
2 and 6, and 2 and 8. These data suggest that BrSIOI and 
BrSIOI work additively in overcoming SI during CO2 treat- 
ment in the C02-sensitive (HA- 11621) fine. No QTL was 
detected at genes known to affect SI stability (MLPK, ARCl, 
or the S-locus), indicating that CO2 sensitivity is determined 
by novel genes in the experimental lines used here. 

Associated gene prediction by in silico comparative 
mapping 

Using the B. rapa genome sequence (Cheng et al., 2011), 
BrSIOI could be mapped to a 569 kb region flanked by InDel 
marker XT05-004 and SSR marker BRMS-034, and BrSIOI 
to a 1469 kb region flanked by SSR markers BRMS-042-2 and 
KBrHl 10I17R. These two regions include 121 and 280 genes 
annotated in the Brassica database (BRAD), respectively 
(Supplementary Tables S3, S4 at JXB online). Comparison 
of the A. thaliana genome with the Brassicaceae genome 
(reviewed by Schranz et al., 2006) suggests that BrSIOI has 
synteny on A. thaliana chromosome 2 and BrSI02 has syn- 
teny on both chromosomes 3 and 4. It is assumed that these 
two QTLs do not have the same genetic origin and could be 
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Fig. 8. QTL analysis results. The solid line indicates the LOD score and the dotted line indicates the QTL threshold (LOD=3.4) 
determined using a 1000 permutation test (P < 0.05). Thex-axis represents each linkage group (cM) and they-axis indicates the QTL 
score. Two QTLs (BrSI02 and 3) are detected in AOS and one in AOS {BrSIOI). Arrows show loci involved in SI stability reported by 
Hatakeyama ef a/. (2010). 
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Table 1 . Summary of CO2 sensitivity QTLs 



QTL 


LG 


Closest marker 


QTL peak (cM)= 


LOD 


R2b 


Additive effect^ 


BrSIOI 


AOS 


>Cr05-004 


83.50 


5.17 


19.30 


0.72 


BrSI02 


AOS 


BRMS-042-2 


60.87 


4.46 


19.00 


0.69 


BrSI03 


AOS 


KBrH110l17R 


41.25 


S.76 


14.50 


0.65 



" QTL peak position, detected by interval mapping, between two markers. 

Amount of phenotypic variation explained by the QTL. 
" Additive effect of the CQg-sensitive HA-1 1 621 allele on RLSICQ2. 



Table 2. Statistical analysis of QTL effect 



QTL 



BrSIOI 



BrSI02 



BrSI03 



Marker" 



KBr B075I01F 



>CT05-004 



BRMS-034 



BRMS-042-2 



BRMS-050 



KBrH110l17R 



BRMS-114 



ss 


3.12 (27f 


3.17 (24)' 




3.24 (20)' 


3.10(23) 




3.09 (22)' 


3.01 (25) 


3.01 (30)' 




SI 


2.34 (53) 


2.36 (51) 




2.42 (52) 


2.43 (60) 




2.47 (63) 


2.45 (62) 


2.30 (54) 




II 


2.02 (26) 


1 .93 (25) 


2.06 (33) 


1 .87 (23) 


1 .84 (22) 


1 .86 (22) 


2.03 (22) 





^ S, CQ2-sensitive HA-1 1 621 allele; /, CO^-insensitive HA-1 1 623 allele. 
Individuals whose genotype was unidentified are excluded. 

Kruskal-Wallis analysis comparing phenotype between genotype groups with individuals in the same groups. 
Significance level: **P <0.01 ; *P <0.05. 



Table 3. QTL association for COg sensitivity 



Group 
no. 


Marker^ 
XT05-CI04 


BRMS-042-2 


No. of 
individuals'' 


IVIean RLSICO2 


1 


SS 


SS 


5 


3.86 






2 


SS 


SI 


13 


3.23 1 








3 


ss 


II 


5 


2.40 








4 


SI 


SS 


13 


3.12 








5 


SI 


SI 


23 


2.37 








6 


SI 


II 


13 


1.63 ■ 








7 


II 


SS 


6 


2.42 






8 


II 


SI 


16 


1.78 






9 


II 


II 


3 


1.71 







" S, CQ2-sensitive HA-1 1 621 allele; /, CO^-insensitive HA-1 1 623 allele. 
Individuals whose genotype was unidentified are excluded. 
Significance level, *P <0.05. 



two independent regions controlling high CO2 sensitivity. 
Based on reciprocal cross results, the CO2 sensitivity trait may 
be controlled by genes expressed in the female organ (Fig. 8). 
A total of 1 2 1 and 280 annotated genes in BrSIOI and BrSI02 
have 103 and 243 homologues in A. thaliana, respectively, and 
54 and 141 of these genes are expressed in A. thaliana pis- 
til (microarray data of carpel at stage 12, http://afFymetrix. 
arabidopsis.info/narrays/search.pl?fl = l&sl=ATGE_37, last 
accessed 14 December 2013, Supplementary Tables S3, S4). 

Genes involved in related biological processes are often 
expressed cooperatively and their co-expression information 
is important for understanding biological systems (Eisen 
et al, 1998). ATTED-II (http://atted.jp/, last accessed 14 
December 2013) is a gene co-expression database useful for 
identifying the potential partners working in the same biolog- 
ical processes (Obayashi et al, 2007). Co-expression analysis 



was performed using ATTED-II with these 195 genes and it 
was found that MAP kinase 6 (At2g43790 in BrSIOI) and 
ethylene overproducer 1 (At3g511770 in BrSIOI) showed the 
strongest co-expression and calmodulin-like 41 (At3g50770 
in BrSIOI) has weak co-expression with cytochrome c oxi- 
dase 10 (At2g44520in BrSIOI) and beta glucosidase 28 
(At2g44460 in BrSIOI). In addition to these co-expressed 
genes, these two regions encode highly homologous family 
member proteins, for example matrixin proteins (At2g45040 
in BrSIOI and At4g 16640 in BrSIOI) and senescence- 
associated proteins (At2g44670 in BrSIOI and At4g 17670 
in BrSI02). All these can be candidate responsible genes, 
although the biological functions of these genes are mostly 
unknown. 



Discussion 

It has been >40 years since Nakanishi et al (1969) first 
reported that SI could be overcome by CO2 and Nakanishi 
and Hinata (1973) demonstrated the applicability of this 
technique to commercial use. Nowadays, seed companies 
have adopted this method to obtain inbred parental seeds of 
crucifer vegetables for large-scale commercial Fi hybrid seed 
production. However, there is still very limited understanding 
of the mechanism by which SI is overcome. 

Lee et al (2001) showed a shrunken and distorted papilla 
cell surface in the C02-sensitive B. rapa cv. Hiratsuka, and 
suggested that these structural changes could cause SI to 
be overcome. The cryo-scanning electron microscopy data 
reported here did not show any structural changes in CO2- 
sensitive or C02-sensitive lines (Fig. 2A). Additionally, pre- 
treatment of non-polhnated pistils with high CO2 gas did not 
cause SI breakdown (Fig. 4). Therefore, a completely differ- 
ent SI breakdown mechanism must be present, at least in 
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the C02-sensitive line. In contrast, massive Ca accumulation 
was observed at the pollen-stigma interface specifically in 
C02-sensitive plants under high CO2 conditions (Fig. 2B). 
Brewbaker and Kwack (1963) were the first to describe 
the need for a high concentration of Ca^"^ for pollen ger- 
mination and pollen tube growth. The high concentration 
of Ca^"^ could be needed for activating pectinase to loosen 
the papilla cell wall, allowing the pollen tube to penetrate 
(Black and Charlwood, 1995), or for keeping the pollen tube 
cell wall rigid enough not to burst (Hepler and Winship, 
2010). Although causal relationships remain unclear, the 
data suggest that CO2 treatment induces a certain compat- 
ible reaction leading to Ca^"^ accumulation at the pollen- 
stigma interface. 

To date, several genetic studies have been performed to 
understand the mechanism of SI breakdown for breeding 
purposes. Niikura and Matsuura (2002) reported that in 
Japanese radish high CO2 sensitivity is controlled by a reces- 
sive gene that governs the construction and/or metabolism of 
the stigma, which reacts to CO2 without any changes in gene 
expression. In contrast, Hyun et al. (2007) reported a domi- 
nant, S-haplotype-linked high CO2 sensitivity phenotype in 
B. rapa. In contrast to these previous reports, Fi plants had 
an intermediate CO2 sensitivity and the F2 population had a 
continuous frequency distribution of RLSICO2 in the present 
study (Fig. 6). These results suggest that in the fines used for 
this study, CO2 sensitivity is a quantitative trait which is con- 
trolled by more than one gene. 

Genetic finkage maps based upon frequency of recom- 
bination in segregating populations are fundamental and 
powerful tools for associating phenotypic trait-specific 
genetic regions. Linkage mapping can be used to understand 
the biological basis of complex traits and to dissect genetic 
determinants underlying the expression of agronomically 
important breeding traits (Paran and Zamir, 2003). Using 
QTL analysis, two major QTLs for high CO2 sensitivity were 
successfully identified (Fig. 8). BrSIOl and BrSI02 had 
similar LOD scores and explained similar amounts of phe- 
notypic variation (19.3% and 19%), and these could be two 
major factors controlling high CO2 sensitivity. Very recently, 
five QTLs associated with stability of SI in B. rapa have been 
identified. Two of them co-locahzed with SLG (A07) and 
MLPK (A03) and the other three were on A02, A06, and 
AlO (Hatakeyama et al., 2010). CO2 sensitivity did not Hnk 
with the S-locus in the present study (Fig. 5) and none of 
the other reported loci co-localized with QTLs detected 
here (Fig. 8), indicating that CO2 sensitivity of the lines in 
this study is determined by novel genes different from those 
known to affect SI stability. Genes in BrSIOl and BrSIOl 
regions have 103 and 243 homologues in A. thaliana, respec- 
tively, and 54 and 141 of these genes are expressed in A. thal- 
iana pistil. In silico comparative analyses identified several 
co-expressing genes and highly homologous genes encoded 
in these two regions. All these can be candidate responsi- 
ble genes; however, to dentify the genes responsible for high 
CO2 sensitivity in the QTL regions in B. rapa more accu- 
rately, it would be necessary to narrow down the regions by 
developing near-isogenic fines (NILs). 



To maintain Fi seed quafity, inbred lines with strong but 
C02-sensitive SI are ideal for Fj hybrid breeding, and it is 
very important to understand the genetic relationships 
between Sl-related genes and CO2 sensitivity phenotypes. 
These results could be useful for the marker-assisted selection 
of parental lines with both stable SI and high CO2 sensitivity. 

Supplementary data 

Supplementary data are available at JXB online. 

Figure S 1 . S-haplotype analysis of F2 plants by PCR-RFLP. 

Figure S2. RLSICO2 in C02-sensitive and C02-insensitive 
lines, and Fi and F2 progeny based on the number of pen- 
etrating pollen tubes after self-pollination under high CO2 
conditions. 

Table SI. Genetic markers and their primers used for link- 
age analysis. 

Table S2. 5'-haplotype segregation in the F2 population. 
Table S3. Annotated genes and Arabidopsis homologues in 
BrSIOl. 

Table S4. Annotated genes and Arabidopsis homologues in 
BrSIOl. 
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